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ABSTRACT 

Let  (x^y^)  , v = be  points  of  interpolation  with 


0 < x. 


We  consider  sequences  {S^}  of 


'1  < * " " < ^ and  let  1 < p <_ 

2ir-per  iodic  functions  which  interpolate  the  data  optimally  in  the  sense  that 


m 


= minimum 


The  main  results,  which  depend  on  the  parity  of  k,  concern  the  asymptotic 

behavior  of  S (x)  as  m tends  to  infinity, 
m 

1.  k = 2n  + 1.  Let  T(x)  be  the  unique  polynomial 


(1) 


U r 

T(x)  = — + l (a  cos  qx  + b sin  qx) 
q=l  q q 


that  interpolates  the  data.  Then 


(2) 


lim  S (x)  = T(x)  uniformly  in  x . 
m 

nr*” 


2.  k = 2n.  Among  all  interpolants  (1)  let  T(x)  be  that  polynomial  such  that 

2 2 

a + b = minimum  . 
n n 

For  this  T (x)  again  (2)  holds. 

AMS  (MOS)  Subject  Classifications:  Primary  41A15;  Secondary  42A12 

Key  Words:  periodic  spline  functions,  trigonometric  polynomials,  interpolation 

Work  Unit  Number  6 (Spline  Functions  and  Approximation  Theory) 

SIGNIFICANCE  AND  EXPLANATION 

We  consider  points  in  the  plane  which  recur  in  a wave-like  pattern.  We 
determine  infinitely  many  curves  passing  through  these  points  and  indexed  by  a 
natural  number  m.  The  index  m indicates  the  smoothness  of  the  curve;  and 
as  m becomes  larger  the  curves  take  on  a limiting  shape  which  is  explicitly 
determined  in  the  report. 
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PERIODIC  INTERPOLATING  SPLINES  AND  THEIR  LIMITS 


A.  S.  Cavaretta,  Jr.  and  D.  J.  Newman 


In  1940  J.  Favard  wrote  "Sur  1 ' interpolation, " a paper  which  has  attracted  some 
attention  during  the  past  few  years  [2],  [3],  [8].  Favard's  work  raises  many  interest- 
ing questions  (e.g.  [3])  and  we  would  now  like  to  address  one  which  appears  in  the  very 
final  paragraph  of  his  paper:  to  wit,  how  does  one  handle  interpolation  of  periodic 
data  by  smooth  m-times  differentiable  functions  and  what  is  the  asymptotic  behavior  of 
these  smooth  interpolating  functions  as  m tends  to  infinity.  The  present  paper  is 
accordingly  divided  into  two  sections.  In  the  first  section  we  approach  the  interpolation 

problem  as  a constrained  minimization  and  obtain  a sequence  S of  spline  interpolants . 

m 

The  asymptotic  behavior  of  this  sequence  is  determined  in  the  second  section. 

We  deal  always  with  uniform  convergence  on  the  period,  while  the  constrained 
minimization  is  posed  in  each  of  the  L^  norms,  1 < p < ®.  However  the  limit  of  the 
sequence  does  not  depend  on  p,  but  rather  is  completely  determined  by  the  given 

data.  The  case  when  p = 2 has  been  settled  earlier  by  M.  V.  Golitschek  [6]  and  also 
by  I.  J.  Schoenberg  [9].  The  new  ideas  necessary  to  handle  the  case  of  p ^ 2,  and 
in  particular  p = ”,  consistute  the  novelty  of  the  present  paper. 
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a □ 


51.  Fix  (x  ,y  ),  v = as  the  points  of  interpolation  and  assume 

v v 

0 1 *1  < X2  < " ‘ < xk  < 2lT  and  k 1 2-  We  denote  by  W^  the  space  of  2n-periodic 
functions  whose  (m-l)th  derivative  is  absolutely  continuous  and  whose  m -derivative 
is  in  the  Lebesgue  space  1,^(0, 2n).  Put 

F = {fe  W^  I f (x  ) = y , v = l,...,k)  . 
p P ' V V 

Our  problem  is  to  determine  solutions  of 

(1.1)  min  ||  f (m)  ||  • 

feF  p 

P 


For  nonperiodic  data  the  results  for  this  interpolation  problem  are  well  known  (4),  f2), 
and  we  merely  adapt  the  effective  approach  of  de  Boor  and  others  to  the  present  case  of 
periodic  data.  When  p = «•  the  result  for  the  perodic  case  reads  as 

Theorem  1.  Problem  (1.1)  has  a solution  S which  is  a periodic  perfect  spline  function 
of  degree  m with  at  most  2 [— ) knots.  More  precisely 

|s(m)  (x)  | = c 


except  at  the  knots  and 


c = min  ||f 
fe  F 


Moreover  if  m > k the  solution  S is 


Remark:  Except  for  the  uniqueness  statement,  this  fact  has  been  observed  by  others, 
notably  B.  D.  Bojanov  [1)  and  W.  Forst  [5] . 

Proof:  We  first  convert  our  interpolation  conditions  into  moment  conditions.  For  every 
integer  j and  v = l,...,k  put 

xv+jk  “ xv  + 2*j 
yv+jk  = yv  • 

Based  on  the  nodes  {x  }°°  , construct  the  B-splines  M of  degree  m - 1 defined 

\)  v=“®  V 

as  the  Peano  kernels  of  divided  differences  of  order  m.  Thus  for  any  f e F^  and  any 
integer  v 


-2- 


ty  • • • • «y  1 

7v  J v+m 


v+m 


Mv(t)f(.ln)  (t)dt 


Now  form  (t)  = J M (t) , v = 1, . . . ,k.  Clearly  the  are  2n-periodic  and 

j=-“  V 3 

linearly  independent. 

Now  for  any  f « f set 


(1.2) 


c 

V 


/ t (t)f (m)  <t)dt, 

A ” 


V = 1. 


,k  . 


We  now  solve  the  constrained  minimization  problem 

2ir 

(1.3)  min|lg||  subject  to  / f (t)g(t)dt  = c , v = l,...,k  . 

" ^ V V 

Let  M denote  the  linear  span  of  . We  view  M as  a k-dimensional  subspace 

of  Lj(0,2n),  and  on  M we  define  a linear  functional  X by  setting 

—Ci  v = 1 , . . . , k . 
v v 

* 

It  follows  from  the  Hahn-Banach  theorem  and  the  identification  of  L,  with  L that 

1 

there  exists  a solution  gt  e L^(0,2-n)  of  -the  above  minimum  problem  (1.3).  Clearly 

II 9*11.  = INI- 

To  determine  the  structure  of  the  function  gt(t),  let  us  suppose  that  X achieves 
its  norm  for  the  function  s e M.  Thus  ||s||^  =1  and  Xs  = ||x||.  Suppose  for  the 
moment  that  the  periodic  spline  s vanishes  on  no  interval.  Then  from  Rolle's  theorem 

Jj 

we  conclude  that  s(t)  has  no  more  than  2 [— 1 zeros  on  its  period.  Observe  also 

(1-4)  IlgJL  - INI  = As  = / s(t)gt(t)dt  < ||gJ|J|s||  = ||gj|„  . 

0 x 

So  equality  must  hold  throughout  and  so 

9.(t)  - ||  g.  ||  sign  s(t)  , 

except  at  2(— ) points  or  less  where  s(t)  = 0. 

Now  if  m >_  k it  follows  from  the  minimum  support  property  of  the  B-splines  M^(t) 

that  s cannot  vanish  on  an  interval.  So  g4(t)  is  in  this  case  uniquely  determined 


by  sign  s(t).  If  m < k and  s(t)  should  actually  vanish  on  an  interval,  we  need 


Thus  the  divided  differences  of  p(t)  vanish;  hence  on  the  nodes  p(t)  must  reduce  to  a 
constant  pQ.  So  S(t)  = S^(t)  - pQ  interpolates  the  data  and  yields  the  desired 
solution  of  Theorem  1. 

For  arbitrary  p,  1 < p < “,  we  obtain  by  the  same  method  as  above  the  following 
proposition.  In  case  p = 2,  the  interpolating  spline  S is  structurally  characterized 
as  the  interpolating  spline  with  knots  at  the  xy,  see  [8];  for  other  p the  structure  of  S 
is  less  transparent,  but  has  been  discussed  to  some  extent  by  Golomb  [7] , also  by 
de  Boor  (3) . 


Proposition:  For  each  p,  1 < p < »,  there  is  a unique  S t Fp  minimizing  ||s(n')||p, 
This  S i£  characterized  by  the  relation 

s(m)  (x)  = Isfx)!'*  1 sign  s(x) 
where  <3  is  conjugate  to  p and  s is  a spline  in  M. 

Before  leaving  questions  of  interpolation,  let  us  recall  some  results  concerning 
interpolation  by  trigonometric  polynomials.  For  k odd,  say  k = 2n  + 1,  there  is  a 
unique  trigonometric  polynomial  T(x)  of  order  n satisfying 

(1,5)  T(xv)  = yv*  v = 1,  ...,2n  + 1 . ‘ * 

It  is  convenient  at  times  to  write 


(1.6) 


an  n 

T(x)  = -j-  + l (a 


<3=1 


cos  qx  + b sin  qx) 


(1.7) 


T(x) 


2n+l 

I y i (x) 

v=l  v v 


where  t^lx)  are  the  Lagrange  functions  for  the  nodes  x ,. 


2n+l  ‘ 


For  k even,  k = 2n,  there  is  a one  parameter  family  of  interpolating  nth  order 
trigonometric  polynomials  with 

(1-8)  T(xv)  = yv,  v = 1 , . . . , 2n 

and  Schoenberg  [8]  has  suggested  singling  out  the  one  with  least  amplitude  for  the  terms 
of  highest  frequency: 

l^flnition:  Among  the  T(x)  olf  the  form  (1.6)  and  satisfying  the  interpolation  condi- 
tions (1.8)  we  determine  the  unique  T(x)  which  satisfies  the  condition  that 

2.2 

a + b = minimum  . 
n n 


This  T is^  called  the  proximal  interpolant. 

2n  x - x 

Consider  the  function  ij>(x)  = ] f sin  


v=l 


Evidently  t|i  vanishes  precisely  at  the 


xv  arx*  so  Tj  (x)  satisfies  the  interpolation  conditions  (1.8)  so  does  T^(x)  + Ai/>(x) 

for  any  value  A.  Using  this  observation  one  can  easily  establish  the  following  lemma  (9] 


Lemma 
of  the 


: The  proximal  interpolant  T(x)  is  uniquely  characterized  as  the  interpolant 
form 


§2.  For  the  data 


(x^, yv> » v = l,...,k,  we  solve  the  problem  in  Theorem  1 for 

every  positive  integer  m.  This  yields  a sequence  IS  (x) } of  periodic  perfect 

m m=l 

splines,  all  of  which  interpolate  the  same  data.  Alternatively,  obtain  as  in  the 
Proposition  a similar  sequence  {S^ (x)  which  minimizes  || S lip'  1 < p < °“.  In  any 

case,  the  following  theorem  describes  the  asymptotic  behavior  of  these  sequences. 

Theorem  2:  Depending  on  whether  k i£  odd  or  even,  let  T(x)  be  the  unique  interpolating 
trigonometric  polynomial  as  determined  in  SI  (1.6)  for  k = 2n  + 1 and  by  (1.9)  for 
k = 2n . Then  for  every  j = 0 , 1 , . . . , we  have 

lim  S<3>  (x)  = T^  (x)  uniformly  in  x . 
m — 

m-w 

52.1.  The  number  of  data  is  odd:  k = 2n  + 1 . 

We  need  to  write  the  Fourier  series  for  our  functions.  Since  we  will  differentiate 
the  resulting  series,  it  is  convenient  to  use  the  complex  form  even  though  all  functions 
concerned  are  real.  Accordingly,  let 


(2.1) 


S (x) 
m 


l 


c 


1SX 

e 

s,m 


= U (x)  + R (x)  ; 
m m 


here  we  have  set 


U (x)  = 
m 

Lemma  2 : For  j = 0,1,2,..., 

Proof : Observe  that 

1 

c = — 
s,m  2n 

so  by  Holder's  inequality 


l c e1EX,  R (x)  = l c e1 
-n  S'm  ro  |s|>n  S'm 


lim  R (x)  = 0 uniformly  in  x. 
m — 

nr*0 


2tt 


2TT 


/ -hi  s'-»(.)  ■ 


-ist 


. _ ) e dt 

. m 27i  ' m 

0 0 


dt 


(-is) 


Since  T also  interpolates  the  given  data,  we  have  from  the  optimality  property  of 
that 


-7- 


<2. 2) 


. <m) 


lp  HIT 


(m) 


Ip  < ""Nip 


where  the  last  inequality  is  that  of  Bernstein-Zygmund.  Hence 

m 

(2.3) 


s,m  — s 


80 

(2.4) 


Clearly  as  m tends  to  infinity,  the  right  hand  side  of  (2.4)  tends  to  0,  and  thus 
the  lemma  is  proved  for  j = 0.  For  arbitrary  j we  find  in  exactly  the  same  way 
(2.5)  l»(i) 


l«.  . I |f 

I s | >n  1 


m-] 


m > j + 1 , 


which  completes  the  proof. 

Using  Lemma  2,  the  case  of  k = 2n  + 1 of  Theorem  2 follows  quite  easily.  Indeed 


and  from 


it  follows 


S (X)  = U (x)  + R (x) 
m m m 


lim  R (x)  = 0 
m 

m-*» 


lim  °m(xv)  = yv'  v = l,...,2n  + 1 
m-*“ 


But  as  yv  = T (x^)  and  each  U^  is  a trigonometric  polynomial  of  order  n,  we  conclude 

lim  U (x)  = T (x) 

m 

or*” 

uniformly  in  x.  Finally  observing  that  the  operation  of  differentiation  is  continuous 
when  restricted  to  trigonometric  polynomials  of  order  n,  we  have 

lim  sj[3\(x)  = lim  U^(x)  + lim  R(^  (x)  = lim  U(^(x)  = (x) 

m-K»  m m-«“  m m-xo  m m 

uniformly  in  x. 

S2  2.  The  number  of  data  is  even:  k = 2n. 

The  even  case  is  more  subtle  and  to  prove  the  result  we  need  the  precise  value  of 


lim  n "s 


(m)  | 


1 < p < 


-8- 


which  is  itself  a result  of  independent  interest.  For  convenience  in  what  follows,  we 
take  the  norm  to  be 


Ilf  II, 


f.  2m  -sp 

[^T  / |f<*>  |Pdxj  1 1 P < 00 


ess  sup  | f (x) | 


P = 


As  a preliminary  observation  in  the  case  k = 2n  we  first  pick  a point  xQ  ^ x^  mod  2n, 
v = l,...,2n.  We  then  prove 

Lemma  3:  Fix  1 < p ® and  let  (S^)  be  the  sequence  of  interpolating  splines  as  ' ‘ 

determined  in  SI.  Then 


(2.6)  S (x.)  = 0(1)  . 

m 0 

Proof:  Recall  that  S (x)  = U (x)  + R (x)  and  from  Lemma  2 (which  still  applies  as  the 

in  iu  m 

parity  of  the  number  of  data  points  plays  no  role) 

lim  R (x)  = 0 . 
m 

m-w» 

Thus  it  is  sufficient  to  show 

(2.7)  n 5 U (x  ) = 0(1)  . 

m mu 

Using  the  Lagrange  interpolation  formula,  we  have 

Um(x)  = j"  (yv  + e^)tv(x)  + T^l0(*) 

where  again  from  Lemma  2 e™  = 0(1),  v = l,...,2n.  Now  for  some  a * 0 and  £, 

fQ(x)  = a sin (nx  + 5)  + lower  terms. 

Then  the  same  argument  which  produced  (2.3)  yields 

2m 

(2.8)  / S (x)sin(nx  + £)dx  = 0(1)  . 

0 ro 

On  the  other  hand 


2x  2x 

/ S (x)sin(nx  + £) dx  = / (U  <x)  + R (x))sin(nx  + £)dx 
0 m 0 

2tt 

= / U (x)sin(nx  + £)dx 
(2.9)  0 m 

2x  2n  2 r 

= / (5|  (y  + em)  i (x))sin(nx  + £)dx  + n / t (x)sin(nx  + £,)dx 

0 v=l  v v v m 0 

= 0(1)  + n ■ air  . 
m 

Together  (2.8)  and  (2.9)  yield  (2.7). 

Taken  in  conjunction  Lemmas  2 and  3 guarantee  that  subsequential  limits  of  the 

sequence  S (x)  exist  uniformly  in  x.  In  particular  from  Lemma  3 we  have  that  the 
m 

trigonometric  polynomials  are  bounded  at  2n  + 1 distinct  points  and  so  the  sequence 

U has  subsequential  limits.  But  from  Lemma  2 it  fellows  that  any  such  limit  (which 
m 

must  be  a trigonometric  polynomial)  is  also  a limit  of  the  corresponding  subsequence 

of  S . We  finish  the  proof  of  Theorem  2 by  showing  that  the  only  such  limiting 
m 

trigonometric  polynomial  is  the  unique  proximal  interpolant  T singled  out  in  §1.  This 
is  proved  by  the  following  two  lemmas,  where  we  argue  for  subsequences  without  indicating 
this  in  the  notation. 

Leim.-.a  4 : Let  { } denote  any  subsequence  and  suppose 

lim  S (x)  = T(x)  = a sin(nx  + u)  + lower  order  terms  . 

. m 

m-**» 

Then 


(2.10) 


lim  inf|Jn  n's*n''(x)||  > -rrr — r rj— 

m~  * P-2llslnxllq 


, 1 1 

l<p<o°  — + — =1 

p q 


Proof:  We  argue  for  m even; 


odd  goes  analogously.  Clearly 


2x  2x  _ 

lim  / S (x)sin(nx  + u))dx  = / T(x)sin(nx  + ui)dx  = xa 
m 0 


m-*»  Q 

On  integrating  by  parts  m times  we  get 

2x 


m 

lim  (-1)2  / 

m t 

nr*«>  n 0 


S(m)  (x)sin(nx  + u>)dx  11  xa  ; 
m 
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r 


and  since 

y-  / S{m)  (x)sin(nx  + w)dx|  < ||  S (x)  ||  ||  sin  (nx  * ui)  ||  “ II s ^ <x>  !L  II  sin  x II  q 

^ m | ~ m p 4 111  F ^ 

the  result  follows. 

Lemma  5 : Suppose  for  some  8 and  T 

U(x)  = B sin(nx  + i)  + V(x) 

satisfies  UU^)  = yv,  v = 1, . . .,2n;  V(x)  a trigonometric  polynomial  of  order  n - 1 


(2.11) 


lim  sup  II  n S 


m"(m)  <x>|| 


p — 2 II  sin  x| 


1 < p < 


m -*  M 11  q 

Proof:  By  a translation  we  can  assume  t = 0.  Because  of  the  minimal  nature  of  S^(x) 
we  need  only  produce  for  any  given  e > 0 an  F(x)  which  interpolates  and  for  which 


If (m)  (x)  ll  < 


p — 2 I isin  x 


f — n (lei  + t)  . 


Also  as  seen  in  §1  1 1 S (x ) II  depends  continuously  on  the  data  and  lienee  the  F we 

' m p 

produce  need  only  interpolate  the  data  (x^,  yy) > v = lf.***2n  within  c and  not 
necessarily  exactly. 

For  q conjugate  to  p,  we  consider  the  function 

iq 


Isin  xl 
sin  x 


A = 


2 ( ||  sin  x ) 


with  its  Fourier  series  given  by 


) c sin  sx  . 
s 

S=1 


Form 


(2.12) 


» s 

F(x)  = l c — U (sx) 
s m 

s=l  s 


for  p near  1-  and  U as  in  the  statement  of  the  lemma . Since 


, 2n  | . iq 
1/  . sin  x , , 

*»  — ] A J ; 1 — sin  x dx  = 1 

1 n sin  x 


we  have 


Th<- 


-11- 


F(x)  = PU  (x)  + 0 — 


Thus  F(x)  satisifes  the  interpolation  conditions  within  an  error  of  0(1  - p ) + 0 
which  is  less  than  e for  p close  to  1-  and  m large,  as  they  will  be  chosen. 

Next  observe  from  the  uniform  convergence  of  (2.12)  when  p < 1,  we  have  for  m 
even  (a  similar  analysis  will  hold  if  m odd) 


m 


(2.13)  F ^ ^ (x)  = £ cgP  u""'  (sx)  = (-l)Z6nm  £ c pS  s*n  nsx  + I c pSV^m'  (sx)  . 

s=l  s=l  s=l  S 

Now  observe  that  the  series  of  the  first  term  is  the  Abel  sum  for 

iq 


7 c sin  snx  = A 

*-  C 


S = 1 


For  the  series  of  the  second  term,  we  have 


I sin  nx 1 
sin  nx 


V(m)  = 0(n  - l)m  and  cg  = 0(1)  , 


and  so  the  series  is 


Thus  (2.13)  is  estimated  by 


0( 


1 - P 


(n  - 1)  ) . 


(2.14) 


F(m)  (X)  = (-nV  A ^-n  — + 0(1  - P)nm  + 0(^—  (n  - l)m)  . 

sin  nx  1 - p 


We  make  the  error  in  (2.14)  less  than  cn  by  first  choosing  p sufficiently  close  to 

1-  and  then  choosing  m sufficiently  large.  Taking  the  L norm  we  then  have 

P 


(2.15) 


But  as 


I n_lnF  (m)  (x)  ||  p £ |6|  A 


+ e . 


l |J_  f2lt  lain  nx  1^ 
1 2 it  l ,p 


" l 'sin  x'q) 


A ||sin  x||P  = — n — : r 

11  "q  2 1|  sin  x j 


0 | sin  nx|*  / ’0  ’ " ” ”"q 

(2.15)  implies  (2.11). 

Now  Theorem  2 follows  immediately  from  Lemmas  4 and  5 and  the  previous  remarks.  The 
proximal  interpolant  T,  by  its  very  definition,  provides  a U for  Lemma  5 with  the 
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smallest  possible  B.  Thus  (2.11)  holds  with  this  B.  So  by  Lemma  4 in  conjunction 

with  (2.11)  no  subsequence  of  S can  converge  to  any  other  trigonometric  polynomial 

m 

other  than  the  proximal  interpolant.  It  follows  that  the  full  sequence  S converges 

m 

uniformly  to  T.  We  then  can  argue  just  as  in  §2.1  that  for  j = 0,1,... 

lim  S^(x)  = T^  (x)  uniformly  in  x . 
m 

m+o' 

We  will  close  with  a corollary  concerning  the  perfect  splines  of  Theorem  1,  the 
case  p = °°. 

Corollary:  Let  T be  the  interpolant  of  our  Theorem  2 and  let  a denote  the  amplitude 

of  the  terms  of  frequency  n in  T.  Then 

, . ~mn  (m)ii  ti 

lim  ns  = - a . 

11  m 1 » 4 

m-*» 

Proof:  This  follows  immediately  from  Lemmas  4 and  5 and  the  observation  that  when 
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